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Abstract. Quantum stabilizer states over Fm can be represented as self-dual 
additive codes over F^2 . These codes can be represented as weighted graphs, 
and orbits of graphs under the generalized local complementation operation 
correspond to equivalence classes of codes. We have previously used this fact 
to classify self-dual additive codes over F4. In this paper we classify self- 
dual additive codes over Fg, Fig, and F25. Assuming that the classical MDS 
conjecture holds, we are able to classify all self-dual additive MDS codes over 
Fg by using an extension technique. We prove that the minimum distance 
of a self-dual additive code is related to the minimum vertex degree in the 
associated graph orbit. Circulant graph codes are introduced, and a computer 
search reveals that this set contains many strong codes. We show that some of 
these codes have highly regular graph representations. 

1. Introduction 

It is well-known that self-orthogonal additive codes over F4 can be used to repre- 
sent a class of quantum error- correcting codes known as binary stabilizer codes [4]. 
Although the binary stabilizer codes have been studied most, several authors have 
considered nonbinary stabilizer codes over finite fields [1] [TTJ [T^ [Tni HH US] , cyclic 
groups [13], and Abelian groups in general [55]. We will focus mainly on codes 
over finite fields, and exploit the fact that a stabilizer code over corresponds 
to a self-orthogonal additive code over . Quantum codes of dimension zero are 
known as stabilizer states, which are entangled quantum states with several possible 
applications. Stabilizer states correspond to self-dual additive codes. It is known 
that such codes can be represented as graphs [TH [55] . It is also known that two self- 
dual additive codes over F4 are equivalent if and only if their corresponding graphs 
are equivalent, up to isomorphism, with respect to a sequence of local complemen- 
tations O m EHl El] . We have previously used this fact to devise a graph-based 
algorithm with which we classified all self-dual additive codes over F4 of length up 
to 12 |7j. Recently, the representation of equivalence classes as graph orbits was 
generalized to self-dual additive codes over any finite field [2] . In this paper we use 
graph-based algorithms to classify all self-dual additive codes over Fg, Fig, and F25 
up to lengths 8, 6, and 6, respectively. We also give upper bounds on the number of 
codes, derived from mass formulas. By using a graph extension technique we find 
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that there are only three non-trivial self-dual additive MDS codes over Fg, assuming 
that the classical MDS conjecture holds. We prove that the minimum distance of 
a self-dual additive code is related to the minimum vertex degree in the associated 
graph orbit. Finally, we perform a search of circulant graph codes, a subclass of the 
self-dual additive codes, which is shown to contain many codes with high minimum 
distance. The highly regular graph structures of some of these codes are described. 



2. Stabilizer states 

Data in a classical computer are typically stored in bits that have values either 
or 1. Similarly, we can envisage a quantum computer where data are stored in 
quantum bits, also known as qubits, i.e., two-level quantum systems. One qubit 
can then be described by a vector |a;) = (^) G C^, where |ap is the probability 
of observing the value when we measure the qubit, and is the probability 
of observing the value 1. More generally, data could be stored in m- level qudits, 
described by vectors from C". Measuring such a qudit would give a result from 
an alphabet with m symbols. In general, this alphabet could be any finite Abelian 
group, but we will focus on the case where the alphabet is a finite field. The m 
vectors \x), x € F^, form an orthonormal basis of C™. 

An error operator that can affect a single qudit is represented by a complex 
unitary mxm matrix, i.e., a matrix U such that UU'' = I, where f means conjugate 
transpose. A state of n qudits is represented by a vector from C™ — C™ (S> ■ ■ • C™. 
Assuming that errors act independently on each qubit, this state is affected by error 
operators described by n-fold tensor products of unitary mxm matrices. In the 
case of qubits (m — 2), we only need to consider errors from the Pauli group, 



X = 



due to the fact that these matrices form a basis of all unitary 2x2 matrices. The 
error X is called a bit- flip, since AT = |a; -I- 1). The error Z is known as a phase- 
flip, since Z \x) — (--1)^ \x). For general qudits that take their values from F^, we 
consider the generalized Pauli group, Vm, also known as the discrete Heisenberg- 
Weyl group. When our alphabet is a finite field, we must have m = p^ , where 
p is a prime and r > 1. The errors contained in the generalized Pauli group are 
shifl: errors, X{a) \x) = \x-\-a), and phase errors, Z{b) \x) — uj^'^rn/p{bx) |^^^ where 



, b S Fm , cij is a complex p-th root of unity, and tr„ 



/p 



F 



F„ is the trace 



function, tr„/p(c) = J2l 
group is generated by 



If m = p is a prime, i.e., r = 1, the generalized Pauli 



Xil) = 



/o 



Vi 



Z(l) = 



0/ 



/I 



Vo 



OJ 



\ 



-7 



where w is a complex p-th root of unity, and I is the identity matrix of specified 
dimension^ The operators X(a) and Z{b) are obtained by taking the a-th and 6-th 
powers of A'(l) and Z{1), respectively. Even if m is not prime, we can still define 



^The set of generators also contains the scalar lu, except for the case m = 2, where it contains 
i, a 4-th root of unity. This overall phase factor can be ignored for our purposes. 
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qudits that take values from the cychc group Zm , and use the same error operators 
as defined above. However, when m is a prime power, we get much better codes 
by using a finite field as our alphabet. When we work with qudits that take values 
from ¥pr, where r > 1, we use the error group {^l^^ Ei \ Ei € Vp} [T6j, i.e., the 
operators are r-fold tensor products of Pauli matrices from the group Vp. The error 
bases that we use are examples of nice error bases |17j . 

Quantum codes are designed to add redundancy in order to protect quantum 
states against errors due to interference from the environment. A code of length n 
and dimension k adds redundancy by encoding k qudits using n qudits. One type 
of code that exploits the fact that the generalized Pauli group forms a basis for all 
possible errors is the stabilizer code |10| . A stabilizer is an Abelian group generated 
by a set oi n ~ k commuting error operators. An error is detected by measuring 
the eigenvalues of these operators. If a state is a valid codeword that has not been 
affected by error, we will observe the eigenvalue +1 for all operators. The quantum 
code, i.e., the set of all valid codewords, is therefore a joint eigenspace of the stabi- 
lizer. If there is a detectable error, some eigenvalues would be different from +1, due 
to the commutativity properties of the generalized Pauli matrices. A stabilizer gen- 
erated by a set of n error operators defines a zero-dimensional quantum code, also 
known as a stabilizer stateio The minimum distance of a zero-dimensional stabilizer 
code is simply the minimum nonzero weight of all error operators in the stabilizer. 
The weight of an error operator is the number of m x to tensor components that 
are different from the identity matrix. A quantum code of length n, dimension fc, 
and minimum distance d, over the alphabet F„i, is denoted an [[n, A:,c?]]m code. 
Stabihzer states are therefore [[n, 0, d]]m codes. If the minimum distance, d, is high, 
the stabilizer state is robust against error, which indicates that it is highly entan- 
gled. Entangled quantum states have many potential applications, for instance in 
cryptographic protocols, or as graph states [23] which can be used as a resource for 
quantum computations. In the next section we will also see that zero-dimensional 
stabilizer codes correspond to an interesting class of classical codes, known as self- 
dual additive codes. 

Example 1. A [[4, 0,3]] 3 stabilizer state is obtained from the stabilizer generated 
by the following error operators. 

X{1) ^X{l)Z{2)(g,I 
X{l)Z{l)(g,X{2) (g)X{l)Z{l)(E)X{l), 
I A(2)Z(2)®X(1)Z(1)®Z(2), 

X(l) (g)Xi2)Z{2)(g>X{2) (g)X{2)Z{2). 

3. Self-dual additive codes 

We can represent a stabilizer state over by an n x 2n matrix {A\ B) [1] . The 
submatrix A represents shift errors, such that = a if X(a) occurs in the j-th 

tensor component of the i-th error operator in the set of generators. Similarly, the 
submatrix B represents phase errors. 



^Stabilizer states could also be called one-dimensional quantum codes, since they are one- 
dimensional Hilbert subspaces. We use the term dimension to mean the number of qudits the code 
can encode. 



4 



Lars Eirik Danielsen 



{A\B) 



Example 2. The matrix corresponding to the stabihzer state in Example [T] is 

/ 1 1 1 2 ^ 
12 1110 10 
2 1 2 1 2 
\122202G2/ 

The matrix [A \ B) generates a code C, and this code is a representation of a 
stabihzer state. The fact that a stabihzer is an Abehan group translates into the 
requirement that C must be self- dual with respect to a symplectic inner product, 
i.e., 

(a I 6) * (a' I 6') = tr„,/p(6 • a' - 5' • a) = 0, V(a | 6), (a' | 6') £ C. 
We define the symplectic weight of a codeword {a \ b) G C as the number of posi- 
tions i where a^, bi, or both are nonzero. (This is the same as the weight of the 
corresponding Pauli error operator.) 

We can also map the linear code of length 2n defined above to an additive code 
over F„j2 of length n. The representation of binary stabilizer codes as self-dual ad- 
ditive codes over F4 was first demonstrated by Calderbank et al. [1] , and generalized 
to qudits by Ashikhmin and Knill [T], and by Ketkar et al. [16j . An additive code, 
C, over Fto2 of length n is defined as an Fm-linear subgroup of F^a. The code C 
contains m" codewords, and can be defined by an n x n generator matrix, C, with 
entries from F„i2 , such that any Fm-linear combination of rows from C is a code- 
wordlfl To get from the stabilizer representation (A \ B) to the generator matrix C, 
we simply take C = A + ljB, where is a primitive element of F^2. The code C 
will be self-dual, C = C^, where the dual is defined with respect to the Hermitian 
trace inner product, — {it 6 F^a \ u * c = for all c G C}. When m — p 
is prime, the Hermitian trace inner product of two vectors over Fp2 of length n, 
u = (ui,U2, ■ • ■ ,Un) and v = (wi,f2, • ■ • ,Vn), is given by 

n 

1=1 

When m = p^ is not a prime, we use a modification of the Hermitian trace inner 
product [16], 



tr 



m/p 



UJ - UJ'"- 

where w is a primitive element of F„2 . 

The Hamming weight of a codeword u, e C, denoted wt(M), is the number of 
nonzero components of u. The Hamming distance between u and v is wt(M — v). 
The minimum distance of the code C is the minimal Hamming distance between 
any two distinct codewords of C. Since C is an additive code, the minimum distance 
is also given by the smallest nonzero weight of any codeword in C A code over F„j2 
with minimum distance d is called an (n,m^,d) code. The weight distribution of 
the code C is the sequence {Aq, Ai, . . . , An), where Ai is the number of codewords 
of weight i. The weight enumerator of C is the polynomial 

n 

Wix,y) = J2A,x"-Y 

i=0 

For an additive code over F„2, all Ai must be divisible by m — 1. 



^For additive codes over F4, each codeword is a sum of rows of the generator matrix. However, 
we also use the name "additive code" in this more general case. 
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Example 3. The stabilizer state in Example [T] corresponds to the following gener- 
ator matrix of a self-dual additive (4, 3^, 3) code. 



C = 



/I l + 2uj 1 \ 

2 l+uj 1 

2 + 2uj l + LU 2uj 

\ 1 2 + 2uo 2 2 + 20)/ 



We define two self-dual additive codes, C and C over F,„2 , to be equivalent if the 
codewords of C can be mapped onto the codewords of C by a map that preserves 
the properties of the code, including self-duality. A permutation of coordinates, 
or columns of a generator matrix, is such a map. Other operations can also be 
applied to the coordinates of C. Let each element a + cjfo G V^i be represented as 
(b) G J^m- We can then premultiply this element by a 2 x 2 matrix. (We could 
equivalently have applied transformations to pairwise columns of the 2n x n matrix 
{A I B).) It was shown by Rains [22 that by applying matrices from the symplectic 
group Sp2(m) to each coordinate, we preserve the properties of the code. (This 
group contains all 2 x 2 matrices with elements in F„j and determinant one.) For 
self-dual additive codes over F4, these symplectic operations can be represented 
more simply as multiplication by nonzero elements from F4 and conjugation of 
coordinates. (Conjugation of elements in ¥p2 maps x to x^.) Combined, there 
are six possible transformations that are equivalent to the six permutations of the 
elements {l,Ci;,cj^} in the coordinate. The corresponding symplectic group is 

where Ai represents multiplication by uj and A2 represents conjugation. Including 
coordinate permutations, there are a total of 6"n! maps for a code of length n. 
For codes over Fg, we observe that Sp2(3) is a group of order 24 generated by 

where Ai represents multiplication by uj'^ and A2 represents the map a + ujb 1-^ 
a + b + Lub. By taking powers of Ai, we see that we are allowed to multiply a 
coordinate by a; € Fg only if xx = 1. However, if we also conjugate the coordinate, 
we may multiply by a; g Fg where xx = 2. Note that conjugation on its own is 
not allowed. The 8 operations just described may be combined with the operations 
represented by A2 and A2 to give a total of 24 operations. In all there are 24"n! 
maps that take a self-dual additive code over Fg to an equivalent code. In general, 
for codes over F^2, the number of maps is | Sp2(m)|"n!. 

A transformation that maps C to itself is called an automorphism of C. All 
automorphisms of C make up an autom,orphism group, denoted Aut(C). The number 
of distinct codes equivalent to a self-dual additive code over ¥^2, C, is then given by 

|Sp,(^ - 



^^^^^^1 . The equivalence class of C contains all codes that are equivalent to C. 
By adding the sizes of all equivalence classes of codes of length n, we find the total 
number of distinct codes of length n, denoted T„. The number T„ is also given by 
a mass formula. The mass formula for self-dual additive codes over F4 was found 
by Hohn [T5| . This result is easily generalized to F^a . 
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Theorem 3.1. 

n t. 



a iAut(c,)i ' 

where tn is the number of equivalence classes of codes of length n, and Cj is a 
representative from each equivalence class. 

Proof. Let M{n, k) be the total number of self-orthogonal (n, m*^) codes. One such 
code, C, can be extended to a self-orthogonal (n, m'^+^) code in m^^""*"'' — 1 ways 
by adding an extra codeword from C^. Each [n^m^'^^) code can be obtained in 
this way from rn^lfe+i) _ ^ different {n,m'') codes. It follows that 

Min,k + l) = M{n,k)-^^^^^. 

Starting with M{n,0) — 1, the recursion gives us the number of self-dual (n, m") 
codes, 

By assuming that all codes of length n have a trivial automorphism group, we get 
the following lower bound on i„, the total number of inequivalent codes. Note that 
when n is large, most codes have a trivial automorphism group, so the tightness 
of the bound increases with n. Also note that this bound is much tighter than a 
bound that was derived from results in graph theory by Bahranigiri and Beigi [2]. 



Theorem 3.2. 

^nLi(^' + i) 

I Sp2(TO)|"n! 

where c = 1 if m is even, and c = 2 if m is odd. 



tr,. > 



Proof. When m is even, the trivial automorphism group includes only the identity 
permutation, and the result follows from Theorem l3.1l When m = is odd, where 
p is a prime, the trivial automorphism group also contains the transformation that 

applies the symplectic operation q ^ p^ ij ^'^ coordinates. This operation 

is equivalent to multiplying each codeword by p — 1, and will therefore map an 
additive code to itself. □ 



It follows from the quantum singleton bound [18\ I22j that any self-dual additive 
code must satisfy 2d < n + 2. A tighter bound for codes over F4 was given by 
Calderbank et al. ^4]. Codes that satisfy the singleton bound with equality are 
known as maximum distance separable (MDS) codes. Self-dual MDS codes must 
have even length, and MDS codes of length two are trivial and exist for all alphabets. 
The only non-trivial MDS code over F4 is the (6,2^,4) Hexacode. Ketkar et al. [16l 
Thm. 63] proved that a self-dual additive {n,m'^,d) MDS code must satisfy n < 
m? + d — 2 < 2m^ — 2. If the famous MDS conjecture holds, then n < m? + 1, or 
n < m? 2 when m is even and d = 4 or d = vi? . Grassl, R5tteler, and Beth [24] 
showed that MDS codes of length n <m+l always exist. 

Self-dual linear codes over F„2 are a subset of the self-dual additive codes. Only 
additive codes that satisfy certain constraints can be linear. Such constraints for 
codes over F4 were described by Van den Nest [5D] and by Glynn et al. [S]. An 
obvious constraint is that all coefficients of the weight enumerator, except Aq, of 
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a linear code must be divisible by — 1, whereas for an additive code they need 
only be divisible by m — 1. 



A graph is a pair G = {V,E) where ^ is a set of vertices and E C V x V 
is a set of edges. Let an m-weighted graph be a triple G = (V, E, W) where W 
is a set of weights from ¥m- Each edge has an associated non-zero weight. (An 
edge with weight zero is the same as a non-edge.) An m-weighted graph with n 
vertices can be represented by an n x n adjacency matrix T, where the element 
^{ij) = W{{i,j}) if {i,j} G E, and = otherwise. We will only consider 

simple undirected graphs whose adjacency matrices are symmetric with all diagonal 
elements being 0. The neighbourhood of v G V, denoted Ny c V, is the set of 
vertices connected to v by an edge. The number of vertices adjacent to v, \Ny\, 
is called the degree of v. The induced subgraph of G on [/ C V contains vertices 
U and all edges from E whose endpoints are both in U. The complement of a 
2- weighted graph G is found by replacing E with V x V — E, i.e., the edges in E 
are changed to non-edges, and the non-edges to edges. Two graphs G = {V, E) and 
G' = {y, E') are isomorphic if and only if there exists a permutation -k oiV such that 
{u, v} £ E <^=> {7r(?i), 7r(w)} G E' . We also require that weights are preserved, i.e., 
^{u,v} = W{Tr{u),T^{v)}- A path is a sequence of vertices, {vi,V2, ■ ■ ■ ,Vi), such that 
{vi, V2}, {v2, V3}, . . . , {vi-i,Vi} G E. A graph is connected if there is a path from 
any vertex to any other vertex in the graph. A complete graph is a graph where all 
pairs of vertices arc connected by an edge. A clique is a complete subgraph. 

Definition 4.1. A graph code is an additive code over Fj„2 that has a generator 
matrix of the form C = T + u)I, where I is the identity matrix, a; is a primitive 
element of F„2, and F is the adjacency matrix of a simple undirected m-weighted 
graph. 

Theorem 4.2. Every self-dual additive code otjer F„2 is equivalent to a graph code. 

Proof. The generator matrix, C, of a self-dual additive code over F„2 corresponds 
to an n X 2n matrix {A \ B) with elements from F^, such that C = A + ujB. We 
must prove that an equivalent code is generated by (F | /), where / is the identity 
matrix and F is the adjacency matrix of a simple undirected m-weighted graph. A 
basis change can be accomplished by {A' \ B') = M{A | B), where M is an n x n 
invertible matrix with elements from F„i. If B has full rank, the solution is simple, 
since B~^{A \ B) = (F' | I). We obtain (F | I) after changing the diagonal elements 
of F' to 0, by appropriate symplectic transformations. Any two rows of (F | I) 
will be orthogonal with respect to the symplectic inner product, which means that 
TI^ — ir"^ = 0, and it follows that F will always be a symmetric matrix. In the 
case where B has rank fc < n, we can perform a basis change to get 



where Bi is a fc x n matrix with full rank, and Ai also has size kxn. Since the row- 
space of {A' I B') defines a self-dual code, and B' contains an all-zero row, it must be 
true that A2BJ = 0. A2 must have full rank, and the row space of Bi must be the 
orthogonal complement of the row space of A2. We assume that Bi = {Bu \ B12) 
where Bu is a k x k invertible matrix. We also write A2 = {A21 \ A22) where 
A22 has size {n — k) x (n — k). Assume that there exists an x e F;^-'= such that 



4. Correspondence to weighted graphs 
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Figure 1: Graph Representation of the (4,3"', 3) Code 



A22X^ = 0. Then the vector v = {0, . . . ,0,x) oi length n satisfies A2v'^ = 0. Since 
the row space of Bi is the orthogonal complement of the row space of A2 , we can 
write V = yBi for some y G Fjj^. We see that yBu = 0, and since Bn has full rank, 
it must therefore be true that y = 0. This means that a; = 0, which proves that 
A22 is an invertible matrix. Two of the symplectic operations that we can apply to 

columns of a generator matrix are ™ " and f \ J V This means that 

y \m — 1 OJ 

we can interchange column i of A' and column i of B' if we also multiply one of the 

columns by m — 1. In this way we swap the i-th columns of A' and B' for k < i < n 

to get {A" I B"). Since Bu and A22 are invertible, B" must also be an invertible 

matrix. We then find B"~'^{A" | B") = (r | /), and set aU diagonal elements of T 

to by symplectic transformations. □ 

Example 4. The matrix from Example [2] can be transformed into the following 
matrix, using the method given in the proof of Theorem 14.21 



(r I /) = 



/Olio 
1001 
1002 
\ 1 2 



1000^ 
0100 
0010 
0001/ 



This means that the stabilizer state from Example [T] is equivalent to the graph code 
generated by C = F + ujI. The graph defined by F is depicted in Fig. [T] 



Note that Theorem 14.21 is a generalization of the same theorem for codes over 
F4 [7], which was proved by Van den Nest et al. [21 . The fact that stabilizer codes 
can be represented by graphs was also shown by Schlingemann and Werner [26^ and 
by Grassl, Klappenecker, and Rotteler [12] . 

We have seen that every m-weighted graph represents a self-dual additive code 
over F,„2 , and that every self-dual additive code over F„j2 can be represented by an 
rn-weighted graph. It follows that we can, without loss of generality, restrict our 
study to codes with generator matrices of the form T + ujI, where F is an adjacency 
matrix of an unlabeled simple undirected m-weighed graph. 

5. Graph equivalence and code equivalence 

Swapping vertex i and vertex j of a graph with adjacency matrix F can be 
accomplished by exchanging column i and column j of F and then exchanging row i 
and row j of F. We call the resulting matrix F'. Exactly the same column and row 
operations map T + ujI to V + ioI, which are generator matrices for equivalent codes. 
It follows that two codes are equivalent if their corresponding graphs are isomorphic. 
However, the symplectic transformations that map a code to an equivalent code do 
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(a) The Graph G (b) The Graph G * 1 



Figure 2: Example of Local Complementation 



not in general produce isomorphic graphs, but we will see that they can be described 
as graph operations. 

It is known that two self-dual additive codes over F4 are equivalent if and only 
if their corresponding graphs are equivalent, up to isomorphism, with respect to 
a sequence of local complementations [3l [9l 1201 121j . We have previously used this 
fact to devise a graph-based algorithm with which we classified all self-dual additive 
codes over F4 of length up to 12 [7 . 

Definition 5.1 ( 3 ). Given a graph G = {V, E) and a vertex v & V , let CV he 
the neighbourhood of v. Local complementation (LC) on v transforms G into G *v 
by replacing the induced subgraph of G on Ny by its complement. 

Theorem 5.2 ( [S] [9l [20l [21] ) . Two self-dual additive codes overWj^, C and C , with 
graph representations G andG' , are equivalent if and only if there is a finite sequence 
of not necessarily distinct vertices (wi, V2, ■ ■ ■ ,Vi), such that G 'i' vi =¥ V2 * ■ ■ ■ * Vi is 
isomorphic to G' . 

The LC operation can be generalized to weighted graphs, and it was first shown 
by Bahramgiri and Beigi [2 that the equivalence of nonbinary stabilizer states over 
Fm, i.e., self-dual additive codes over F„2, can be described in terms of graph 
operations 

Definition 5.3 (f5]). Given an m-weighted graph G = {V,E,W) and a vertex 

V d V, weight shifting on w by a G F^ \ {0} transforms G into Go^ v by multiplying 
the weight of each edge incident on w by a. 

Definition 5.4 (^). Given an m-weighted graph G — {V,E,W) and a vertex 

V V, generalized local complementation on v by a G F™ \ {0} transforms G into 
G *a V. Let F and F' be the adjacency matrices of G and G *a v, respectively. Then 

Theorem 5.5 ([2]). Two self-dual additive codes overWjn^, C and C' , with graph 
representations G and G' , are equivalent if and only if we get a graph isomorphic 
to G' by applying some finite sequence of weight shifts and generalized local comple- 
mentations to G. 



^Bahramgiri and Beigi [2] only state their theorem for Fm where m is prime, but the result 
holds for any finite field, as their proof does not depend on m being prime. 
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(a) The Graph G (b) The Graph G 1 



Figure 3: Example of Weight Shifting 




(a) The Graph G (b) The Graph G *a 1 



Figure 4: Example of Generalized Local Complementation 



A proof of Theorem 15.51 was given by by Bahramgiri and Beigi |^ , as a gener- 
alization of the proof given by Van den Nest et al. [SJ for self-dual additive codes 
over F4. 

Definition 5.6. The LC orbit of a weighted graph G is the set of all non-isomorphic 
graphs that can be obtained by performing any sequence of weight shifts and gen- 
eralized LC operations on G. 

Theorem 5.7. The minimum distance of a self-dual additive {n, m", d) code is equal 
to 6 + 1, where d is the minimum vertex degree over all graphs in the associated LC 
orbit. 

Proof. A vertex with degree d — 1 in the LC orbit corresponds to a codeword of 
weight d, and we will now show that such a vertex always exists. Choose any graph 
representation of the code and let G = (F | /) be the corresponding generator 
matrix. Find a codeword c of weight d generated by G. Let the i-th row of G be 
one of the rows that c is linearly dependent on. Apply symplectic transformations 
to the coordinates of the code such that c is mapped to c' with 1 in coordinate 
n -\- i, and with in all other of the last n coordinates. Since we do not care about 
changes in the corresponding first n coordinates, as long as the symplectic weight 
of c is preserved, there will always be transformations that achieve this. Apply the 
same transformations to the columns of G, and then replace the i-th row with c', 
to get G'. Note that the right half of G' still has full rank, so we can transform G' 
into a matrix of the form (F' | /) by Gaussian elimination, where the symplectic 
weight of the i-th row is d. Finally, we set all diagonal elements of F' to zero 
by appropriate symplectic transformations. Vertex i of the graph with adjacency 
matrix F' has degree d — 1. □ 
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6. Classification 

It follows from Theorem 1 5 . 51 that two self-dual additive codes over F„2 are equiv- 
alent if and only if their graph representations are in the same LC orbit. The LC 
orbit of a graph can easily be generated by a recursive algorithm. We have used the 
program nauty |19j to check for graph isomorphism. 

Let Gn,m be the set consisting of all non-isomorphic simple undirected connected 
rn-weighted graphs on n vertices. Note that connected graphs correspond to inde- 
composable codes. A code is decomposable if it can be written as the direct sum of 
two smaller codes. For example, let C be an {n,Tn^,d) code and C an (n',m" ,d') 
code. The direct sum, C®C' — {u\\v \ u G C,v ^ C'}, where || means concatenation, 
is an {n + n',m"^" , minjd, d'}) code. It follows that all decomposable codes of 
length n can be classified easily once all indecomposable codes of length less than 
n are known. 

The set of all distinct LC orbits of connected m-weighted graphs on n vertices is 
a partitioning of Gn,m into in,m disjoint sets. in,m is also the number of indecom- 
posable self-dual additive codes over ¥jy^2 of length n, up to equivalence. Let Ln,m 
be a set containing one representative from each LC orbit of connected m-weighted 
graphs on n vertices. The simplest algorithm for finding such sets of representatives 
is to start with the set Gn,m and generate LC orbits of its members until we have a 
partitioning of Gn,m- The following more efficient technique is based on a method 
described by Glynn et al. [9 . Let the m" — 1 extensions of an m-weighted graph 
on n vertices be formed by adding a new vertex and joining it to all possible com- 
binations of at least one of the old vertices, using all possible combinations of edge 
weights. The set En,m, containing i„_i^,„(m"~^ — 1) graphs, is formed by making 
all possible extensions of all graphs in Ln^i^m- 

Theorem 6.1. i„_m C En,m, 'i-e-, the set En^m will contain at least one represen- 
tative from each LC orbit of connected m-weighted graphs on n vertices. 

Proof. Let G = (V, E, W) e Gn,m, and choose any subset U C V oi n — 1 vertices. 
By doing weight shifts and generalized LC operations on vertices in U, we can 
transform the induced subgraph of G on [/ into one of the graphs in Ln~i.m that 
were extended when En.m was constructed. It follows that for all G G Gn.m, some 
graph in the LC orbit of G must be part of En^m- CH 



The set En^n will be much smaller than Gn,m, so it will be more efficient to search 
for a set of LC orbit representatives within En^m- Another fact that simplifies our 
classification algorithm is that weight shifting and generalized local complementa- 
tion commute. This means that to generate the LC orbit of a weighted graph, we 
may first generate the orbit with respect to generalized local complementation only, 
and then apply weight shifting to the resulting set of graphs. 

Using the described techniques, we were able to classify all self-dual additive 
codes over Fg, Fig, and F25 up to lengths 8, 6, and 6, respectively. Table [T] gives the 
values ofin^m, the number of distinct LC orbits of connected m-weighted graphs on n 
vertices, which is also the number of inequivalent indecomposable self-dual additive 
codes over F^2 of length n. The total number of inequivalent codes of length n, 
tn, is shown in Table [2] together with lower bounds derived from Thcorem l3.2l The 
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Table 1: Number (i„.m) of Indecomposable Codes of Length n over F„2 



n 



«n,2 



^n,3 ^n,4 ^n,5 



1 1 

2 1 

3 1 

4 2 

5 4 

6 11 

7 26 

8 101 

9 440 

10 3,132 

11 40,457 



1 
1 
1 
3 

5 

21 
73 
659 



1 
1 
1 
3 
6 

25 



1 
1 
1 
3 
7 

38 



12 1,274,068 



numbers i„ are easily derived from the numbers in by using the Euler transform [27j , 



Tables [3l HI and [5] list by minimum distance the numbers of indecomposable codes 
over Fg, Fig, and F25. A database containing one representative from each equiv- 
alence class is available at http://www.ii.uib.no/~larsed/nonbinary/. For the 
classification of self-dual additive codes over F4, we refer to previous work [7^, and 
the web page http://www.ii.uib.n0/~larsed./vncorbits/ 

Note that applying the graph extension technique described previously is equiv- 
alent to lengthening [IJ a self-dual additive code. Given an (n, m", d) code, we add 
a row and column to its generator matrix to obtain an (n-f 1, m""*"^, d!) code, where 
d' < d+1. If follows that given a classification of all codes of length n and minimum 
distance d, we can classify all codes of length n + 1 and minimum distance d -I- 1. 
All length 8 codes over Fg have been classified as described above. By extending 
the 77 (8,3^,4) codes, we found 4 (9,3^,5) codes, and from those we obtained a 
single (10, 3""^°, 6) code. Assuming that the MDS conjecture holds, there are no self- 
dual additive MDS codes over Fg with length above 10. This would mean that the 
three MDS codes with parameters (4, 3'^,3), (6,3^,4), and (10, 3^°, 6) are the only 
non-trivial self-dual additive MDS codes over Fg. The (6,3^,4) and (10, 3^°, 6) are 
constructed as circulant codes in Section [71 A generator matrix for the (4, 3^, 3) 
code is given in Example [H In fact, a (4, m^, 3) code, for any m > 3, is generated 




d\n 




t 



"n 



by 



/ w 1 1 \ 

\ Q bj a 
\ 1 a uj J 
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Table 2: Total Number {tn,m) of Codes of Length n over F„2 



n 


tn,2 


tn,3 




^n,5 


1 


1 


1 


1 


1 


2 


2 


2 


2 


2 


3 


3 


3 


3 


3 


4 


6 


7 


7 


7 


5 


11 


13 


14 


15 


6 


26 


39 


44 


58 


7 


59 


121 


? 


? 


8 


182 


817 


> 


> 21,161 


9 


675 


> 9,646 


> 458,993 


> 38,267,4t 


10 


3,990 


> 2,373,100 






11 


45,144 








12 


1,323,363 








13 


> 72,573,549 









Table 3: Number of Indecomposable Codes of Length n and Distance d over Fg 



d\n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


2 


1 


1 


2 


4 


15 


51 


388 


? 


? 


3 






1 


1 


5 


20 


194 


? 


? 


4 










1 


2 


77 


? 




5 
















4 


? 


6 


















1 


All 


1 


1 


3 


5 


21 


73 


659 


? 


? 



Table 4: Number of Indecomposable Codes of Length n and Distance d over Fie 



d\n 2 3 4 5 6 

2 1 1 2 4 16 

3 12 6 

4 3 
5 

6 

All 1 1 3 6 25 



where a G F^ \ {0, 1}. This code has weight enumerator W{l,y) = 1 + 4(m^ — 
I)y3 + (m2-3)(m2- 1)2/4. 

There are four (9, 3^, 5) codes, all with weight enumerator W{1, y) = 1 + 252y^ + 
1176y^ + 3672y'^ + 7794y^ + 6788y^. None of these are equivalent to circulant graph 
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Table 5: Number of Indecomposable Codes of Length n and Distance d over F25 



d\n 2 3 4 5 6 

2 1 1 2 4 21 

3 1 3 11 

4 6 
5 

6 

All 1 1 3 7 38 



codes. The generator matrices are: 





LU 


2 


2 


1 











2 


2 ^ 




/ 


UJ 


2 


2 


2 


2 


2 


2 


2 


1 \ 




2 




1 











1 


2 


1 






2 
















2 


1 


2 




2 


1 


UJ 








2 


2 


1 


1 






2 










2 





1 





2 




1 










2 


1 








1 






2 








LU 


1 





1 


2 
















2 


LU 


2 


1 


1 


1 






2 





2 


1 


UJ 


1 


2 


2 


1 










2 


1 


2 




1 


1 









2 











1 


LU 





1 


1 







1 


2 





1 


1 


LU 





2 






2 


2 


1 


1 


2 





a; 


2 


1 




2 


2 


1 





1 


1 





LU 


2 






2 


1 





2 


2 


1 


2 




1 


V 


2 


1 


1 


1 


1 





2 


2 






I 


1 


2 


2 





1 


1 


1 


1 


LU 1 


/ 


LU 


2 











2 


1 


2 


1 \ 






a; 





2 





2 


2 





2 


2 \ 




2 


to 





2 


1 


2 





1 














2 


2 








2 





1 










UJ 





1 


2 





1 


1 






2 


2 










1 





2 










2 





LU 





1 





1 


1 









2 





UJ 








1 


2 


1 







1 


1 





LU 


1 


2 


2 


1 






2 











LU 





2 


1 


1 




2 


2 


2 


1 


1 


LU 


1 


2 


1 






2 





1 








LU 


1 


1 


2 




1 











2 


1 


LU 


1 


1 









2 





1 


2 


1 




1 







2 


1 


1 


1 


2 


2 


1 


LU 


1 






2 





2 


2 


1 


1 


1 


LU 


1 


\ 


1 





1 


1 


1 


1 


1 


1 


UJ J 




I 


2 


1 





1 


1 


2 





1 


UJ 1 



The two (7,3'', 4) codes, the three (6,4^,4), and five of the six (6,5^,4) codes are 
equivalent to circulant graph codes generated in Section [71 The last (6,5^,4) code 
has weight enumerator W{\,y) ~ l + 360y'*+ 3024?/^ + 122402/^ and generator matrix 

/ w 2 2 1 \ 

2 cj 1 4 1 

2 0^011 

1 1 w 3 4 ■ 

4 1 3 w 3 

\ 1 1 4 3 w / 

7. Circulant graph codes 

It is clearly infeasible to study all self-dual additive codes of lengths much higher 
than those classified in the previous section. We therefore restrict our search space 
to the m codes over F„2 of length n corresponding to graphs with circulant 

adjacency matrices. A matrix is circulant if the i-th row is equal to the first row, 
cyclically shifted i — 1 times to the right. We have performed an exhaustive search 
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Table 6: Highest Found Minimum Distance of Codes over F^2 of Length n 



n\m 


2 


3 


4 


5 


2 


2 


2 


2 


2 


3 


2 


2 


2 


2 


4 


2 


3* 


3* 


3* 


5 


3 


3 


3 


3 


6 


4 


4 


4 


4 


7 


3 


4 


4 


4 


8 


4 


4 


4 


4 


9 


4 


5" 


5 


5 


10 


4 


6 


6 


6 


11 


5^ 


5 


6 


6 


12 


6 


6 


6 


6 


13 


5 


6 


6 


7 


14 


6 


6 


7 


8 


15 


6 


6 


7 


7 


16 


6 


6 


8 


8 


17 


7 


7 


8 


9 


18 


8* 


8 


8 


10 


19 


7 


8 






20 


8 


8 






21 


8* 


8 






22 


8 


9 






23 


8 


9 






24 


8 


9 






25 


8 








26 


8 








27 


9^ 








28 


10 








29 


11 








30 


12 









of such graphs, the result of which is summarized in Table [6l This table shows the 
highest found minimum distance of self-dual additive codes over various alphabets. 
A code with the given minimum distance has been found in our search, except for the 
cases marked *, where a better code is obtained in some other way and does not have 
a circulant graph representation! and cases marked s, which are not circulant, but 
obtained by a trivial shortening [8] of a longer circulant code. Minimum distances 
printed in bold font are optimal according to the quantum singleton bound. If n is 
even and the quantum singleton bound is satisfied with equality, we have an MDS 
code. 

We here give the first row of a circulant generator matrix for those codes classified 
in Section [6] that are equivalent to circulant graph codes. There is a unique (6, 3^, 4) 
code with weight enumerator W{l,y) = 1 + 120j/'* + 240y^ + 368?/^ generated by 



^See the web page http : //www . codetables ■ de/ 1 for details on how codes over F4 of length 18 
and 21 can be obtained. 
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(wOlllO). There are two inequivalent (7, 3'', 4) codes generated by (cjllOOll) and 
(w022220), both with weight enumerator W{1, y) = l + 70y''+3362/5+812t/*^ + 968y^. 
There is a unique (10, 3^", 6) code with weight emunerator W(l, y) = 1 + 1680y^ + 
2880y^ + 14040y*^ + 22160^^ + 18288yi° generated by (w012111210). There are 
three inequivalent (6,4^,4) codes with weight enumerator W{l,y) = 1 + 225y^ + 
1080J/5 + 2790/ generated by (wOlllO), (wOlalO), and (wOla^lO), where a = is 
a primitive element of F4. There are five inequivalent (6, 5^, 4) codes generated by 
(wOlllO), (a;01210), (w02220), (wl0201), and (wl2221), all with weight enumerator 
W{1, y) = l+ 360y^ + 3024y5 + 12240/. 

For circulant graph codes of higher length that are optimal according to the 
quantum singleton bound, we find that all codes of the same length have the same 
weight enumerator. In the list below, we give the first row of one generator matrix 
for each weight enumerator. 

• (7,4^4), (wllaall), 

W{1, y) = l + 105/ + 1008y^ + 4830/ + 10440/. 

• (9,49,5), (wOOlaalOO), 

W{1, y) = 1 + il^y'' + 3780/ + 23220?/^ + 88155/ + 146610/. 

• (10,410,6), (cjOlOalaOlO), 

W{1, y) = l + 3150y*^ + 18000/ + 111375/ + 366000/ + 550050yi°. 

• (11, 4", 6), (wOOallllaOO), 

W(l, / = 1 + 1386/ + 13860/ + 99495/ + 505560/ + 1511598yi° 
+ 2062404J/". 

• (7,5^4), (wOllUO), 

W{1, y) = l + 140/ + 2184/ + 17080/ + 58720?/^. 

• (9,5^5), (0)00211200), 

W(l, y) = l + 504?;^ + SAQOy^ + 84240y^ + 507420?/^ + 1352560t/9. 

• (10,51", 6), (^001222100), 

W{1, y) = l + 5040y^ + 54720y^ + 508680?/^ + 2704560j/9 + 6492624t/i°. 

• (11, 5", 6), (^0012222100), 

W(l, y) = l + 1848/ + 31680/ + 370260/ + 2977480/ 
+ 14282664t/i° + 31164192?/". 

• (13,513,7), (wOlOOllllOOlO), 

W(l, / = 1 + 6864/ + 118404/ + 1538680/ + 148674242/1° 
+ 972223202/11 + 388930 7762/1^ + 7180186562/1^ 

• (14,51^^,8), (^1011331331101), 

W{1, 2/) = 1 + 72072/ + 816816/ + 104744647/1° + 90 6 7 9 6 802/" 
+ 5445367202/^^ + 20104418882/1^ + 34464939842/". 
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• (17,51^9), (wOOlOlllOOlllOlOO), 

W{1, y)^l + 97240z/ + 1633632^^° + 24504480?/" + 296652720y^2 
+ 27336204002;" + 18749403360yi'' + 89994568992y^5 
+ 269984494620?/i^ + 381154477680?/i^. 

• (18,51*, 10), (wl2134242124243121), 

W{1, y) = 1 + 1050192yi° + 11456640y" + 180442080?/^^ 

+ 1964813760yi^ + 16877613600?/" + 107991522432?/^^ 

+ 485972877960yis + 1372155934320y" + 1829541554640^^^ 

As mentioned in the introduction, stabilizer codes can be defined over any Abelian 
group, not only finite fields. For comparison, we also generated circulant codes over 
Z4. As expected, the minimum distance of these codes are much worse than for 
codes over Fig. We found a (7, 4^, 4)-code over Z|, but for all other lengths up to 
16, the best minimum distance was equal to the best minimum distance of codes 
over F4 of the same length. 

Gulliver and Kim J13] performed a computer search of circulant self-dual additive 
codes over F4 of length up to 30. Their search was not restricted to graph codes, 
so our search space is a subset of theirs. It is interesting to note that for every 
length, the highest minimum distance found was the same in both searches. This 
suggests that the circulant graph code construction can produce codes as strong as 
the more general circulant code construction. Besides a smaller search space, the 
special form of the generator matrix of a graph code makes it easier to find the 
minimum distance, since any codeword obtained as a linear combination of i rows 
of the generator matrix must have weight at least i. If, for example, we want to 
determine whether a code has minimum distance at least d, we only need to consider 
combinations of d or fewer rows of its generator matrix. 

Circulant graphs must be regular, i.e., all vertices must have the same number of 
neighbours. We have previously discovered O [6] that many strong circulant self- 
dual additive codes over F4 can be represented as highly structured nested clique 
graphs. Some of these graphs are shown in Fig. \E\ For instance. Fig. I5bl shows a 
graph representation of the (12,2^^,6) "Dodecacode" consisting of three 4-cliques. 
The remaining edges form a Hamiltonian cycle, i.e., a cycle that visits every vertex 
of the graph exactly once. Notice that all graphs shown in Fig. [5] have minimum 
regular vertex degree, i.e., each vertex has d—1 neighbours, where d is the minimum 
distance of the corresponding code. 

We have discovered some new highly structured weighted graph representations 
of self-dual additive codes over Fg and Fig. Fig. [6] shows two interconnected 5- 
cliques where all edges have weight one, and a 10-cycle where all edges have weight 
two. The sum of these two graphs, such that no edges overlap, corresponds to the 
(10, 3^°, 6) code. Up to isomorphism, there is only one way to add a Hamiltonian 
cycle of weight two edges to the double 5-clique, since there cannot be both weight 
one and weight two edges between the same pair of vertices. The first row of a 
circulant generator matrix corresponding to this graph is (a;012111210). 

As a second example. Fig. [7] shows two pairs of 4-cliques, each of which is con- 
nected by a length 8 cycle, and two 16-cycles where all edges have weight a and 
a^, respectively, where a = is a primitive element of F4. The (16,4^^,8) code 
generated by (wOa^lalOOOlOOOlala^) corresponds to a sum of these three graphs. 




Figure 6: Two Graphs Whose Sum Corresponds to the (10, 3^°, 6) Code 




Note that the vertices of the graphs corresponding to circulant (10, S"'^", 6) and 
(16,4^^,8) graph codes have degree higher than d — 1. We have tried to obtain 
similar graph representations for other codes in Table[6l but without success. Many 
of the circulant graph codes have vertex degree much higher than d — 1, for in- 
stance the (14, 51-*, 8) code generated by (^1221202021221), and the (18,5l^l0) 
code generated by (wl2134242124243121). 
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